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Abstract 

Many quantum field theoretical models possess non-trivial solutions which are stable for topolog- 
ical reasons. We construct a self-consistent example for a self-interacting scalar field-the quantum 

m 

(or dressed) kink-using a two particle irreducible effective action in the Hartree approximation. 
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This new solution includes quantum fluctuations determined self-consistently and nonperturba- 
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tively at the 1-loop resummed level and allowed to backreact on the classical mean-field profile. 
This dressed kink is static under the familiar Hartree equations for the time evolution of quantum 
fields. Because the quantum fluctuation spectrum is lower lying in the presence of the defect, 



Qh! the quantum kink has a lower rest energy than its classical counterpart. However its energy is 



higher than well-known strict 1-loop results, where backreaction and fluctuation self-interactions 
are omitted. We also show that the quantum kink exists at finite temperature and that its profile 
broadens as temperature is increased until it eventually disappears. 
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Topological defects are perhaps the best known nonlinear solutions of quantum field 
theories. They are explicitly spatially inhomogeneous yet localized and have finite energy 
configurations (sometimes involving defect-antidefect pairs) Moreover topological 

defects, particularly extended ones (strings, domain walls), have configurational entropies 
that increase exponentially with their size. For these reasons defects are important vehicles 
of disorder. They can lead to phase transitions not predictable in their absence, such as 
Kosterlitz-Thouless |3], or become the mesoscopic traces of high energy disorder, as hap- 
pens in defect formation under sudden quenches. Defects can also lead to nonperturbative 
screening effects, e.g. of static magnetic fields by 't-Hooft-Polyakov monopoles 0. 

Because defects are localized solutions carrying quantized topological numbers they can 
often be thought of as generalized quantum particles. In fact in certain circumstances topo- 
logical solutions dominate the phenomenology and should then be thought as the fundamen- 
tal excitations of a new (often more tractable) model, related to the original by a duality 
transformation. Such changes of perspective have been the source of valuable insights into 
the nonperturbative properties of models of condensed matter, high-energy particle physics 
and quantum gravity. 

While the existence and stability of topological solutions follows from general considera- 
tions about the symmetries and dimensionality p| of a given model, the detailed properties 
of topological solutions, e.g. their size and energy, depend on the back-reaction from the 
background of quantum (and thermal) fluctuations in which they live. For example the 
difference in energy between the dressed defect and its classical counterpart will affect esti- 
mates for the probability of its spontaneous nucleation via quantum or thermal fluctuations 
The same rest energy shift also affects the defect's dynamical response. Only a self- 
consistent treatment, where mass corrections occur to all orders in the loop expansion h\ 
can begin to probe the motion at non-zero momentum and long times 2}. In some cases 
fluctuations may even stabilize nontopological classical field profiles Q], but it is not known 
whether 1-loop results survive a self-consistent treatment. 

Incorporating the effects of fluctuations is technically quite difficult, which explains in part 
why most previous work has relied on a strict loop expansion, neglecting either backreaction 
on the classical mean-field profile or self-coupling of the fluctuation degrees of freedom. 
For example, calculations carried out to 2- loop order appear in Ref. and a variational 
self-consistent calculation without backreaction on the mean field is performed in Ref. [l^l • 



Since the corresponding "solutions" are not self-consiste nt, in particular they will not be 
static under non-equilibrium quantum field evolutions |U|. Fermionic backreaction has 
been considered in the calculation of the sphaleron energy, which is relevant for estimating 



baryon number violation 12, 



13 1 . This computation showed the sphaleron energy to be 



quite insensitive to such corrections, however quartic self- couplings of the scalar field were 
neglected Q]. It is precisely these nonlinear field selfinteractions which we wish to consider 
here. 

To start addressing the above issues we construct here an example of a self-consistent 
quantum topological defect. For simplicity, and for the benefit of analytical cross-checks 
on our numerical results at both the classical and 1-loop quantum levels, we analyze the 
familiar 4 kink in 1+1 dimensions. There has also been a lattice Monte Carlo study of 
the quantum kink [14], providing a benchmark against which we can measure the efficacy 
of the self-consistent Hartree approximation. The methods developed here promptly gen- 
eralize to other cases. Self-consistent spherically symmetric droplets, in particular, will be 
discussed elsewhere We include quantum effects at the 1-loop resummed level with 

backreaction on the mean field by applying the familiar Hartree approximation to a IP I 
effective action [16]. This technique places the two-point correlation function on the same 
level as the one-point function, or mean field, thus backreaction is included as naturally as 
the nonperturbative resummation of fluctuations. At this level of approximation, we obtain 
a local, nonlinear eigenvalue-type equation for the quantum fluctuations. The next level of 
improvement, at 2-loop resummed, would lead to a nonlocal integro-differential problem, 
which would require methods of solution different from the ones developed below. Specifi- 
cally we take the \<f) 4 model with potential 

4 y a 

Then the classical kink solution has the well known analytical form 

0(x) = -^tanh(-^r). (2) 

To set up the quantum problem we separate, as usual, the quantum field into a mean 
if = ((p) and a fluctuation operator-valued field tj) ((0) = 0), such that <f>(x,t) = <p(x,t) + 
ip(x,t). The Hartree (or 1-loop resummed) equations for the mean field and the connected 



(Wightman) two-point function G(x,t; x' ,t') = (tp(x,t)ip(x' ,t')) then become 

[□ - At 2 + \<p 2 {x, t) + 3XG{x, t)] <p(x, t) = 0, (3) 
[D + X (x,t)}G(x,t;x',t') = 0, (4) 
x (x, t) = -/j, 2 + 3\cp 2 (x, t) + 3XG(x, t), (5) 

with □ = d 2 — d 2 and G(x,t) = G(x = x',t = if). We seek a static solution to equations 
(jniEJ), obeying boundary conditions <p(x — > +oo) = — ip(x — > — oo) = ip , where <p is one of 
the degenerate vacuum expectation values of the quantum field. These boundary conditions 
imply that the mean field must cross zero at some xq, which defines the center of the kink. 

Eqs. dUEI) can be converted into a (non-linear) eigenvalue problem by the familiar proce- 
dure of decomposing G(x, t; x', if) in a complete orthonormal basis of mode fields, which we 
shall denote by {i/))~(x,t)}. The procedure has been detailed extensively elsewhere and 
we do not repeat it here. The result is 

G(x, t; x', = i( N ( k ) + VM*, t)P k W, + N(k)r k (x, t)M*', 0} > (6) 
k 

where N(k) is occupation number in this basis, which becomes the Bose-Einstein distribution 
ns(k), if the system is in thermal equilibrium N(k) — > ns{k) = (e hulk / T — l) , with ojk the 
energy eigenvalue associated with the eigenvector ipk(x, t). In particular, the equal-point 
function acquires the very simple form 

G(x, t) = J2 [2iV(*0 + 1] M*, tWkfa t). (7) 
k 

Next we specify a basis {ip^x, £)}, together with the scalar product 

[i/)k(x, t).i/> q {x, t))=i / dx {ip* k {x, t)d t ip q {x, t) - ip g (x, t)d t ip* k {x, t)} . (8) 



Because we seek a static solution the self-consistent x( x ) is time independent. Then each 
ipk{x,t) can be separated into the product of functions of time and space as 



Mx,t) = J^e-^ t g k (x), (9) 
and the inner product (jHJ) can be written in terms of the gkix) as 

{9k{x).g q {x)) = dx g* k (x)g q (x), (10) 



which is familiar from quantum mechanics. With these definitions orthogonality in {g^} is 
equivalent to orthogonality in {ipk}- Moreover the completeness relation 

^9k{x)g* k {x') = 5 D (x-x), (11) 

k 

must be satisfied as a result of the canonical equal-time commutation relations obeyed by 
the original quantum field. Finally the explicit self-consistent eigenvalue problem reads 

[-Vl + X(x)]g k (x)=<4g k (x), (12) 

with the eigenmodes obeying periodic boundary conditions in the volume x G {—L,L}: 
gk(L) = gk(—L), and solved in conjunction with equation (j3J) for (p. In practice we discretize 
the fields on a spatial lattice of size N and volume L. This provides a regularization of the 
quantum field theory. In addition we have to devise a renormalization scheme to render all 
results finite as the lattice spacing a = L/N is taken to zero, i.e. in the continuum limit. 

The 2-point function G, as defined in Eq. (J7J), contains an ultraviolet logarithmic diver- 
gence in the continuum limit. We specify the fluctuation physical mass and field expectation 
value over the trivial vacuum. These choices implicitly define the bare mass parameter fi 2 
as a function of the ultraviolet cutoff A = it /a so as to absorb this divergence. Explicitly 

Xo = 2m 2 , Lfl = m 2 /\; — ► fi 2 = m 2 + 3XG , (13) 
f A dk h 

where G = , (14) 

J -a 2tt 2y/k 2 + 2m 2 

Substitution of fi 2 , hence defined, into expression (JHJ) for x( x ) now leads to finite results, 
independent of A, as the continuum (and infinite volume limit) are taken. 

The only remaining divergences occur in the calculation of the energy. The total energy 
can be split into a classical piece E c \ and the quantum (and thermal) fluctuation contribution 
E q . Taking expectation values we can write them as 

Eq = [dxJ2 {Ms) h V ' + X(x)] [1 + 2iV(*0] - \^G 2 {x) 

k 

= 5m - J dx G 2 {x) + J2 UJ nN(k n ), (15) 



Er\ = dx 



(V ^ ))2 + V-/^7A) 2 



(16) 
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The analysis of divergences is familiar at 1-loop, where it has been resolved in at least two 
distinct ways: by direct analysis of the divergences in a large momentum expansion 0, 17 1 



and via the normal ordering of the fluctuation Hamiltonian operator, which is sufficient in 
1 + 1 d,_ to extract the pn yS1 ea! Mt e cc—c. fi fl Q. Because tne y a r e 
implemented very differently we have employed both methods as a check. 

The most straightforward procedure, in our opinion, follows from direct analysis of the 
divergences in a large momentum expansion. For consistency we should start by substituting 
fi 2 defined in (JT3j) into (jl5H16|) . which eliminates all logarithmic divergences, up to an infinite 
constant. This constant results from the trace 5m = ^2 n w n /2. For large enough n, oj n — ► 
cu® = ^/k 2 , + 2m 2 ~ \k n \, resulting in a leading quadratic divergence. This divergence is also 
characteristic of the trivial vacuum energy and should be cancelled against it. With these 
prescriptions we can write the renormalized E q as 

K = ^EK( 1 + 2A ^)) - u n] - \g I dx ix(x) - xo] -\^jdx [G(x) - G } 2 , (17) 

n J J 

The classical energy piece E c \ is now to be evaluated as in (|lfi|). with the substitution 
/i 2 — > m 2 . 

Alternatively a finite energy expression for the trace in 5m can be constructed by suit- 
ably manipulating the fluctuation Hamiltonian operator. The normal ordered Hamiltonian 
is written in terms of a complete set of creation and annihilation operators in the kink 
background, which are then transformed via a Bogoliubov transformation to the trivial vac- 
uum basis, where vacuum expectation values are computed. The residual contribution is 
finite and results from the existence of a non-zero particle number as seen back in the kink 
background mode basis. The final result is 



5m = — Tr 
4 



(0-0o 



(18) 



O 

where O = [— <9 2 + x( x )} an d Oq = O with x( x ) = Xo, corresponding to the trivial vacuum. 
The advantage of this procedure is that it is independent of the exact form of xi x ) i n 0fe 
and is thus equally valid at strict 1-loop order as well as in the resummed self-consistent 
approximation. In practice we cast this expression in the explicit form 

,2 



4 

hi 



(19) 



There remains a final ambiguity in expressions ()17l - 119|) . namely the prescription of the 
upper (ultraviolet) cutoff in the mode sums. In the presence of a non-trivial background 



one expects a spectrum that contains a finite and discrete set of bound modes, lying below 
a continuum of scattering states. Physically it is the advent of these bound states that 
produces the leading (in magnitude) contribution to the energy shift of the quantum kink 
relative to the classical solution. Even for the lattice regularized theory, where the number of 
modes is finite and their eigenvalues discrete, there are at least two seemingly natural ways 
to choose the upper cutoff: including contributions up to a given mode number (say N), or 
alternatively, setting a cutoff in the magnitude of the eigenvalues. These two procedures 

essence of this problem was analyzed 
in detail by Rebhan and van Nieuwenhuizen j2l|, who showed that only a mode number 
cutoff leads to consistent results. By including all of the N modes, bound and unbound, in 
our numerical treatment of these sums, we effectively implement this prescription without 
difficulty 

We are now left with the task of solving the interacting self-consistent problem numer- 
ically. We use a standard relaxation routine J22I to solve Eq. (jH)), together with a set of 

r~i 

LAPACK packages [23( to solve the real symmetric eigenvalue problem Eq.(j3]|3J). At each 
step the fluctuation eigenproblem is solved to produce a new intermediate Gi nt (x), which 
we combine with the old value to produce a new trial G n ew = (1 — l)G \d + jG int , until 
convergence in (p and G(x) is reached, up to a specified precision. The adjustable parameter 
< 7 < 1 controls the size of the update and can be adapted to optimize convergence. 

Both mean-field and fluctuation update steps can be problematic and require some further 
discussion. The system retains an overall translational mode and iterative procedures where 
updates of if and G are staggered can excite it, thus degrading convergence. To eliminate 
this problem we pin the kink's center xq : f(xo) = to the mid-point of our spatial lattice. 
The full profile for the mean field is then obtained through the symmetry ip(— (x — xq)) = 
~if(x—xo). To strict 1-loop order there is a well-known fluctuation zero mode associated with 
the eigenvector go{x) ~ d x <p, regardless of the specific mean-field background. This follows 
immediately by taking a spatial derivative through the (classical) Eq. Q and comparing the 
result to the 1-loop eigenvalue problem. This construction ceases to apply once fluctuations 
are made self-consistent, as their potential is shifted by a term ~ 3XG(x). Neglecting 
back-reaction on <p(x) and assuming this effect to be small, knowledge of G(x) allows us to 
estimate the energy shift of the former zero mode as —>■ 3A J dx gl(x)G(x)go(x), where 
go(x) is the 1-loop zero mode eigenvector. This compares well with with our numerical 



approximation Energy cJq u> 2 w| 



classical 



0.9428 



1-loop 



0.4716 








2 



2m 2 



1-loop resummed 0.7184 0.50m 2 1.71m 2 2m 2 



Hartree 



0.6634 0.33m 2 1.29m 2 1.90m 2 



TABLE I: Numerical results for the kink rest energy in the classical, strict 1-loop, 1-loop resummed 
excluding and including (Hartree) back-reaction of quantum fluctuations on the mean field tp. Here 
we have taken N = 800, L = 20 and A/m 2 = H = 1. We also show the three lowest lying eigenvalues 
for each of these approximations. 

results, even if A is not small. It predicts uo 2 = 0.513m 2 vs. uo 2 , = 0.504m 2 , with the 
parameters of Table HI shown below. Thus we see that no fluctuation zero-mode persists in 
the self-consistent problem. We must nevertheless emphasize that the dressed kink retains 
an overall translational invariance at no energy cost (as the energy functional is independent 
of where the kink is placed) but that this transformation requires a simultaneous translation 
of the mean-field profile and its self-consistent fluctuations. 

The advantage of studying the kink is that several analytical results are known, which 
we can use as checks on our numerical methods. The classical kink energy takes the simple 
form E = ^p^y-- With our choice of parameters this is E = 0.9428, which coincides with 
the energy of our numerical solution, as shown in Table HJ 

At 1-loop order the spectrum of fluctuations about the classical kink is also known exactly. 
The spectrum consists of two discrete modes with eigenvalues Uq = 0, uj\ = |m 2 lying below 
a continuum with io\ = k 2 + 2m 2 . The corresponding (unnormalized) eigenvectors are 



with z = mx/y/2. We verified that to this order our numerical solutions coincide with these 
results. Fig. shows the eigenvectors obtained numerically for low lying modes. Corre- 
sponding eigenvalues and kink rest energies computed at different levels of approximation 
are summarized in Table HI The 1-loop results for the mass shift 5m computed via two 
different methods (jTHj) and (JT7J) agree well with each other (—0.4712 vs. —0.4742) and also 
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FIG. 1: Four of the lower eigenmodes (n = 0, 1, 2, 6, clockwise) in the 1-loop approximation (red), 
1-loop resummed in G (green) and full problem including backreaction on ip (black). 
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with the known analytical result at the 1-loop order, 5m = —0.4711, see Refs. 
Beyond this point no exact analytical results are known. 

As a first step beyond 1-loop, we may consider solving the self-consistent problem Eg. (|T2*Jl 
without backreaction on ip. This case has been studied by Boyanovsky et al. |ipj who 
devised a variational approach, based on the knowledge of the 1-loop perturbative spectrum. 
They assume that a solution for the fluctuation modes can be constructed by treating the 
inverse width of the kink potential A and the asymptotic fluctuation mass M as variational 
parameters and using the eigenvalue equation to determine a relationship between them. If 
we impose the constraints appropriate to our choices, A = m/y/2 and M = y2m, we can 
use the solutions of Ref. to predict the form of the renormalized G(x) to be 



G(x) 



coshW with F(,) = -tan-^), 



(22) 



where rj = Aj \f M 2 — A 2 and the lowest eigenvalue to be to 2 , = M 2 — A 2 . With our choice of 
parameters this would give 7] = 0.408, F(jj) = 0.025 and Uq = 3/2m 2 . Instead our numerical 
solution gives an amplitude for G of 0.226, and uJq = 0.50m 2 , as we have discussed above. 
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FIG. 2: The profiles of the quantum kink, its classical counterpart and G(x) for A = h = 1 

The reason for the discrepancy is that the self-consistent mode profiles, although reminiscent 
of their 1-loop forms, display a different width both from the classical kink profile and from 
each other. This allows them to considerably lower their associated eigenvalues and leads 
to a much larger G than in Ref. |icj]. In other words, the self-consistency in G creates a 
repulsive effect (a decrease in the depth of the attractive potential provided by the classical 
kink) in ()12j) and therefore leads to larger bound-state eigenvalues, which are the principal 
source of the energy shifts. As a consequence the quantum kink energy is now higher than 
at 1-loop, but still smaller than its classical value. Results for the first few modes are 
summarized in Table 

The fully self consistent solution, including backreaction on the mean-field, is shown in 
Fig. |21 The most important qualitative change to the fluctuation spectrum is that the 
lower lying modes become more tightly bound relative to the 1-loop resummed case, see 
Table HI and that a third (shallow) bound state appears. The kink mean-field profile is 
now broader than its classical counterpart. Although the self-consistent mean- field solution 
is not perfectly fit by a tanh(xA), we can quantify this effect by producing a best value 
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FIG. 3: Quantum kink profiles at different temperatures. The kink mean field profile broadens 
with temperature and eventually disappears at T ~ 0.65m. 

A = m/2.07 [c.f. A = mj\[2 for the classical kink]. The broadening of the mean-field 
raises the classical energy but this energy cost is more than compensated for by the lowering 
of the fluctuation spectrum. Allowing the mean-field profile to relax to a more favorable 
configuration leads to the lowering of the self-consistent kink energy. Nevertheless the effects 
of the fluctuation self-interactions still result in a total energy higher than at 1-loop. The 
response of the mean-field and fluctuations taken together, self-consistently, can be thought 
of as a screening effect of the defect by the fluctuations. The mass of the self-consistent kink 
in this approximation is in excellent agreement with the lattice results of Ref. [l^ . 

Finally we investigate the effects of turning on temperature, see Fig. 01 Numerically this 
is best achieved by slowly increasing T and successively constructing solutions at higher 
temperature from previous results at lower T. Fig. |3] shows that the kink broadens as the 
temperature is increased. At T ~ 0.65m the solution starts probing a new local minimum 
around (p = 0, which is an inadequacy of the approximation (the Hartree approximation 
predicts a first-order thermal phase transition for the <p A model, which is known to have 



a second order transition). However, as the kink energy now becomes comparable to the 
temperature, it is plausible that a stable kink no longer exists; our numerical methods 
invariably converge to the trivial ip = solution. In a model with a true second order 
phase transition one may expect self-consistent defects to persist and broaden all the way 
up to T = T c , where they may melt away (their core size diverges) and cease to be localized 
objects. 

In summary we constructed a fully self-consistent solution for a topological defect, in- 
cluding the back-reaction of quantum fluctuations at the level of the Hartree approximation. 
We have shown that at zero temperature the self-consistent dressed kink is lighter than its 
classical counterpart-the kink is still attractive and leads to a lower lying spectrum of fluc- 
tuations relative to the vacuum-but heavier than the well known 1-loop result, due to the 
repulsion among self-consistent fluctuations. At finite temperature the effects of fluctuations 
are enhanced by a thermal population of low lying modes and the kink becomes ever broader 
as the temperature is increased until it vanishes. 
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